Alternative criterion for two-dimensional wrapping percolation 
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Based on the differences between a spanning cluster and a wrapping cluster, an alternative cri- 
terion for testing wrapping percolation is provided for two-dimensional lattices. By following the 
Newman-Ziff method, the finite size scaling of estimates for percolation thresholds are given. The 
results are consistent with those from Machta's method. 
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I. INTRODUCTION 

In a square lattice each site is independently either 
occupied with probability p, or empty with probability 
1 — p. A cluster is a group of occupied (nearest) neigh- 
bor sites pj. With more and more sites being occupied, 
clusters grow larger and larger. For the square lattice 
with free boundaries, once a cluster grows large enough 
to touch the two opposite boundaries, the spanning per- 
colation has occurred. This cluster is called as a spanning 
cluster. While for the lattice with periodic boundary con- 
ditions, it is somewhat more difficult to detect a wrapping 
cluster, which wraps around the lattice. 

Since the very early work on three dimensional poly- 
mers [1[ , percolation has found a variety of uses in many 
fields. Typically in the study of networks people 
keep casting great interests. Recent studies of percola- 
tion in physics mainly involve various functional materi- 
als or components, such as optical lattices [7|], magnetic 
materials or ordinary materials [8rJl0j , and nanocompos- 
ites [lT| - tT3j . Very recently, the study on explosive per- 
colation arose "explosive" interests |L4l - [l7| . Many other 
interesting works on percolation theory and its wide ap- 
plications are collected in some books [l|, [HI, [l9j] . 

In all these studies of percolation, a naive question is 
that how to tell the onset of percolation. The percola- 
tion threshold (or critical probability) p c is the value of 
p for which a spanning cluster (in the case of free bound- 
aries) or a wrapping cluster (periodic boundary condi- 
tions) appears for the first time [2(| Hl| . For a periodic 
lattice, once a wrapping cluster appears, say along the x- 
direction, if we cut the lattice along y-direction to let the 
x-boundaries be open, there exists (at least) one spanning 
cluster. In turn, when a spanning cluster first appears 
on a free boundaries lattice, if we connect the two oppo- 
site open boundaries to let it be a periodic lattice, the 
wrapping cluster does not necessarily appear. Therefore, 
wrapping percolation always happens later than span- 
ning percolation for the lattices with the same linear di- 
mension L but with different boundary conditions. The 
difference between a wrapping cluster and a spanning 
cluster gives us a clue to build an alternative criterion 
for wrapping percolation. Once a spanning cluster oc- 
curs, we amalgamate those clusters (excluding the span- 
ning cluster) touching the boundaries, build connected 
relations between the spanning cluster and the amalga- 



mated clusters. If finally, two ends of the spanning cluster 
connect to each other via an amalgamated cluster, then 
a wrapping cluster appears. The data processing after 
obtaining a wrapping cluster follows the Newman-Ziff al- 
gorithm 20] exactly. 

The prevalent criterion for wrapping percolation intro- 
duced by Machta et al. (22[ has been described in detail 
for bond percolation in Newman-Ziff algorithm. It can be 
easily extend to site percolation, once the displacements 
of neighboring sites to the same root site differ by an 
amount other than zero or one lattice spacing, the clus- 
ter wrapping has occurred. Taking site percolation as an 
example, the core idea of Newman-Ziff algorithm is that, 
starting with an empty lattice, a percolation state can be 
realized simply by adding sites one by one to the lattice, 
a sample state with n + 1 occupied sites is achieved by 
adding one extra randomly chosen site to a sample state 
with n sites. An important technique in the Newman- 
Ziff algorithm is the application of binomial distribution. 
Taking the number of occupied sites n as "energy" , if we 
can get a set of measurements {Q n } in the microcanon- 
ical ensemble, then the observable Q(p) in the canonical 
ensemble is given by 
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Q(p) = £ )p n (l-p) N - n Q, 
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Such observable Q can be probability of cluster wrap- 
ping, mean cluster size, and correlation length etc. The 
main advantage of this technique is that, the continuous 
observable Q(p) for all p can be determined from discrete 
N + 1 values of Q n . 

For convenience, we use the same symbols as in 
Newman-Ziff algorithm. There are four types of probabil- 
ity Rl{p) of cluster wrapping on periodic square lattice 
of L x L sites, Iv^ is the probability of cluster wrap- 
ping along either the horizontal or vertical directions, or 
both; R^ , around one specified axis but not the other 
axis; R^\ in both horizontal and vertical directions; R^f 1 

(v) 

and R T ' around the horizontal and vertical directions, 



respectively. For square lattice, R 
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from which we get the values of R^' and R L h ^ by measur- 
ing only the values of R L e ^ and R^ , or vice versa. Given 
the exact value of Roo{Pc), the solution p of the equation 



Rl(p) = Roo(Pc 



(4) 



gives a very good estimator for p c . The solution p con- 
verges to p c according to 



P-Pc 



-11/4 



(5) 



In this way we get the corresponding values of p c from 
Roo{Pc), Raoipc) and Rx (p c ) respectively. Since the 
wrapping probability R^\p) is nonmonotonic, we use 
the position of its maximum to estimate p c , instead of 

the value of R^(p c )- 

The CPU time Tl is related to the statistical errors a Pc 

according to a Pc ~ T L 1 ^ 2 i 1 / 4 . For finite L and a specific 
algorithm, Tl depends only on n, which represents the 
number of runs of the algorithm, instead of the number 
of occupied sites. To keep the same statistical errors on 
systems of different size, we take appropriate n to fulfill 

t l ~Vl. 

In Sec. |ll] we describe our algorithm. In Sec. IIIII we 
give the values of p c from the algorithm, and compare it 
with those from Machta's method. In Sec. IIVI we give 
our conclusions. 



II. THE ALGORITHM 

Since our criterion for wrapping percolation originates 
from the difference in wrapping cluster and spanning 
cluster, the central task of the method is how to han- 
dle those clusters touching the boundaries. Each run of 
the algorithm starts from a lattice with periodic bound- 
ary conditions. In the following, we take the site wrap- 
ping percolation along x-direction as an example. The 
whole algorithm can be separated into two parts, the 
spanning process and the wrapping process. In the span- 
ning process, we check the states of a newly occupied 
site and its (nearest) neighbour sites previously occupied. 
If one site and its neighbour happen to be on the left 
boundary and right boundary respectively, we call them 
as a pair of occupied "quasi-neighbour" sites. Along x- 
direction, if we shift the periodic boundary conditions to 
the open boundary conditions, two sites of any pair of 
quasi-neighbour sites are no longer neighboring to each 
other. In Fig. 1(a) six pairs of quasi-neighbour sites on 



the square lattice 6 x 6 are (0,5), (6,11), (12,17), (18,23), 
(24,29), (30,35). Two occupied quasi- neighbor sites be- 
long to their respective clusters. Therefore, pairs of oc- 
cupied quasi-neighbour sites correspond to pairs of their 
respective clusters. No implementations are required to 
these pairs of clusters or these occupied quasi-neighbour 
sites. For non-quasi-neighbour sites, if the two sites point 
to the same root site (belong to the same cluster), we 
need do nothing; otherwise, if the two sites belong to 



different clusters, we must merge them into a single clus- 
ter. We do exactly what we us ually do [23| before the 
appearance of a spanning cluster [29[ on a square lattice. 
In other words, the periodic boundary conditions along 
x-direction are suppressed temporarily in the spanning 
process. 

Once a spanning cluster occurs, we turn to the wrap- 
ping process immediately, which aims to find a wrapping 
cluster. Obviously, this process is the core of the present 
algorithm. A spanning cluster has at least two ends 
(left end and right end), which touch the left boundary 
and the right boundary respectively. When the periodic 
boundary conditions along x-direction are recovered, by 
across the boundaries, the two ends of the spanning clus- 
ter could be connected, which implies the appearance of 
a wrapping cluster. This process can be finished in three 
steps. 

In the first step, we check the pairs of clusters touching 
the boundaries by scanning the occupied quasi-neighbour 
sites row-by-row. If two clusters that a pair of occupied 
quasi-neighbour sites respectively belong to own different 
root sites and neither of them is the spanning cluster, we 
merge them into a single cluster. The root site of any of 
them can be chosen as the root site of the merged cluster. 
After the implementations on all the non-spanning clus- 
ters touching the boundaries, what we do next is to build 
relations between the spanning cluster and the (merged) 
non-spanning clusters touching the boundaries. We may 
meet three interesting cases when we scan pairs of clus- 
ters touching the boundaries once more. The simplest 
case is that the two clusters are the same spanning cluster 
(see Fig. 1(a) as an example). We simply add a pointer 
from the spanning cluster to itself. The second case is 
that, a non-spanning cluster does not possess a pointer 
pointing to the spanning cluster, we could add a pointer 
from the former to the latter, and label the boundary 
touched by the non-spanning cluster. For the third case, 
a non-spanning cluster already has a pointer pointing to 
the spanning cluster and previously touches the opposite 
boundary, we add a pointer from the spanning cluster 
to the non-spanning cluster. Two examples for the third 

respectively. 



case are showed in Figs. 1(b) and 1(c) 



For 

other cases, we need do nothing further. In the final 
step, a simple while statement is used to check whether 
the spanning cluster points to itself or not by following 
the pointers added in the second step. 

Thus the wrapping process of our algorithm can be 
summarized as follows. 

(1) Amalgamate pairs of non-spanning touching- 

boundary clusters. 

(2) Add pointers from the amalgamated non-spanning 

touching-boundary clusters to the spanning cluster, 
or in turn, add a pointer from the spanning cluster 
to a non-spanning touching-boundary cluster. 

(3) By following a succession of pointers added above, 

check if we can get from the spanning-cluster to 
itself. 
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FIG. 1: Examples of wrapping clusters with occupied sites 
shaded. Pairs of clusters (excluding the spanning cluster) 
are merged at the first step of scanning the occupied quasi- 
neighbour sites, e.g. (18, 23) in (b) and (30, 35) in (c). In 
(a), the occupied quasi-neighbour sites (12, 17) belong to the 
same cluster. In (b) or (c), two ends of the spanning cluster 
are connected via a merged cluster. 



If no wrapping cluster appears, we occupy one addi- 
tional site on the previous lattice and repeat the span- 
ning process and the wrapping process. Since we aim 
to check the reliability of our algorithm via the calcula- 
tion of percolation threshold only, we halt the algorithm 
once the percolation along both x and y directions are 
detected. The program for our algorithm written in C is 
available online [2J|. 



III. THE RESULTS 

All computations are implemented on a desktop PC 
with CPU clock speed 2.6 GHz and memory 1.96 GB. 
We fix n, runs of the algorithm, to 2 x 10 6 for the lattice 
with L = 256. The computation takes CPU time about 
37 hours. The other values of n for L equal to 128, 64 
and 32 are respectively chosen to ensure T L ~ y/L. The 
total n is about 2.7 x 10 8 . The finite siz e scal ing of p c 
for square lattice L x L is showed in Fig. 2(a) and the 
values of p c obtained from four different probabilities Rl 
are listed in Table |U 
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TABLE I: The values of p c for infinite lattice. Here, P and M 
represents the present method and Machta's method respec- 



tively; p c ' , p, 
wrapping probabilities R 
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To check the reliability of our method, and as a com- 
parison, a similar computation within the framework of 
Machta's method is implemented. For L = 256, n is also 
taken as 2 x 10 6 , CPU time is about 4.5 hours, which is 
only one eighth of that in our method. In other words, 
the computation time in units of hours is Tm — 0.28\/~L 
for Machta's method, while for our method Tp = 2.2>\/Tj. 
The other values of n corresponding to L are chosen by 
reference to this short CPU time, and the total n in 
Machta's method is about 8 x 10 7 . In comparison with 
the simulation of more than 7 x 10 9 separate samples in 



FIG. 2: The values of p c on square lattices of L x L obtained 
from the probabilities of cluster wrapping along one axis but 
not the other (solid squares), both axes (solid circles), one 
axis (solid upward-pointing triangles), and either axis (solid 
downward-pointing triangles), (a) Our method; (b) Machta's 
method. 



the work of Newman-Ziff, the statistical errors in this 
work are larger. The finite size scaling and the results of 
p c are respectively showed in Fig. 2(b) and Table fl] 



Besides the difference in computation time, there are 
tiny differences in the results. Obviously, from Table HI 
one can see that the data of p c and p c e ^ are better 



(smaller errors) than the data of p c b ^ and p c n) in our 
method, while in the Machta's method, on the contrary, 
the data of p c b ^ and p c h ^ are better. This difference pos- 
sibly arises from the larger values of pi^ (and therefore 
larger values of p c b ^ and p c h ^) at small L in our method, 
which can be seen in Fig.O Except these differences, our 
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results coincide with that from Machta's method quite 
well. 



IV. THE CONCLUSIONS 

A wrapping cluster appears, if two ends (touching two 
opposite boundaries respectively) of a spanning cluster 
are connected to each other by across the boundaries. 
The results of percolation threshold in our method are 
as good as that from Machta's method, and are consis- 
tent with the published estimates of the square site per- 
colation threshold j25| - {27| . In comparison with Machta's 
method, although our method is not competitive in com- 
putation time, it provides a unified method for wrapping 
percolation and spanning percolation. This work helps us 
to choose an appropriate method for the further study of 
some kind of spatially correlated percolation model [28[ , 
where the percolation thresholds have not yet been well 
determined. 

The lengthy computation time in our method mainly 
comes from two aspects of our specific algorithm. In 
the spanning process, the clusters are amalgamated by 
suppressing the periodic boundary conditions along x and 
y directions respectively; after a spanning cluster appears 



(either along x or y direction), a wrapping process for 
testing a wrapping cluster is implemented. While in the 
Machta's method, a wrapping cluster could be tested in 
the process of merging clusters on a square lattice with 
full periodic boundary conditions. 

Either although the values of percolation threshold 
could be obtained from spanning percolation on a lat- 
tice with open boundary conditions or from wrapping 
percolation on a lattice with periodic boundary condi- 
tions, different boundary effects are definitely covered in 
the results. With our method, it is possible to give a 
direct computation of different boundary effects. 

The present algorithm differs from the Newman-Ziff 
algorithm only in the criterion for wrapping percolation. 
Without question, the former can be extended to the 
calculation of p c for three-dimensional percolation on the 
cubic lattice in the same way as that of the latter. 
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